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Abstract 

The local Hi mass function (HiMF), like the optical luminosity function, is an impor- 
tant observational input into models of cosmology and galaxy evolution. It is a helpful 
framework for assessing the number density of gas rich dwarf galaxies, which are easily 
missed in optically selected galaxy samples, as well as for determining the cosmological 
density of neutral gas at the present epoch. For H i masses larger than 10^ Mq the HiMF 
is determined with reasonable accuracy and the same function is obtained from both op- 
tical and Hi selected samples of galaxies. However, the faint tail below Mhi < 10'' Mq is 
still ill-constrained and leaves room for a population of gas rich dwarfs or free floating H i 
clouds which hypothetically could contribute significantly to the local gas density. De- 
termining the faint tail far below Hi masses of lO"^ Mq will be a great challenge for the 
future. 



1 Introduction 

The Hi mass function (HiMF), the Hi equivalent to the better known optical luminosity func- 
tion (LF) defines the number of galaxies per cubic Mpc as a function of Hi mass Mhi. This 
function is an important constraint in theories of cosmology and galaxy evolution. The shape of 
the local HiMF tells us how the neutral gas in the nearby Universe is distributed over galaxies 
of different masses. For example, semi-analytical models of galaxy formation require detailed 
measurements of the HiMF and the LF to test the predictions (e.g., [0]). In contrast to the 
LF, the HiMF should also constrain the number density of gas rich dwarfs and low surface 
brightness (LSB) galaxies. It has been shown that these objects are easily missed in optical 
surveys which are used to evaluate the LF ||, ^], but they could contribute a significant part of 
the HiMF. Also the often speculated upon, but yet unidentified class of free floating Hi clouds 
are definitively missed in LFs but are hypothetical contributors to the HiMF. 



Besides the shape of the HiMF, also the normahzation is of great importance to cosmology. 
The integral over the function yields the cosmological neutral gas density of the local Universe. 
This measurement anchors the estimates of the neutral gas densities at higher redshift deter- 
mined from damped Lya systems, seen in absorption in the spectra of background quasars. 
These systems are generally considered to be the high z equivalents of the cold gas disks of 
present day spiral galaxies. The exact measurement of the gas density as a function of redshift 
brackets theories of structure formation and also helps in the understanding of the physical 
processes of gas cooling and star formation. An important question is whether the LSB galax- 
ies and the gas rich dwarfs make a significant addition to the cosmological neutral gas density 
at the present epoch. 



2 How To Measure an HiMF 

Different methods can be employed for evaluating an HiMF. The largest variety lies in the 
sample of galaxies that is used to derive the function. Basically two classes of galaxy samples 
can be distinguished: optically selected galaxy samples and H I selected galaxy samples. This 
section briefly describes both methods. 

2.1 HiMF Based on Optically Selected Galaxies 

An HiMF can easily be derived from optical luminosity functions when a conversion factor 
from optical luminosity to H I mass is used. Suppose the LF can be described by a Schechter 
function of the form 



= <!>*[-) exp-(^-jd(^-j, (1) 

where L* is the characteristic luminosity that defines the knee in the LF, 0* is the normalization, 
and a is the value that defines the slope of the faint end. By using a relation between optical 
magnitude M and H I mass, for example log Mm = A — BM, an HiMF can be derived: 

The low-mass slope of the HiMF is related to the faint-end slope of the LF as am = («opt + 
1)^ — 1. Rao & Briggs [^] applied this method by using different LFs and conversion factors 
for different morphological types and added up the results to obtain a global Mm distribution 
function. The HiMF that they derived is well-fit by a Schechter function with a ~ —1.1. 
Intuitively one would believe that a trend of increasing H I richness for low luminosity galaxies 
{B < 0.4) implies a steep low mass end of the mass function. However, Rao & Briggs already 
stressed that for realistic values of B (~ 0.35 fl^l)) Q^hi ~ Q^opt, especially for Oopt in the range 
-1.0 to -1.3. 

In Fig. m an example of this method is shown. The luminosity function for various morpho- 



logical types taken from Marzke et al. [23| is shown in the left panel. The right panel shows 
the derived Hi mass functions, using different values of A and B for the different morpholo- 
gies. It is clear that the HiMF is dominated by spiral galaxies, although at the very low mass 
end the irregulars seem to overtake the number densities. However, this dominance occurs at 
luminosities below Mb = —14, the magnitude limit of the Marzke et al. sample. 

A slightly modified form of this method of finding an HiMF for optically selected galaxies is 
to use the measured H I fiuxes for a sample of galaxies with well-understood selection criteria. 




Figure 1: Left panel: Luminosity function for various morphological types taken from Marzke 
et al. [Q. Right panel: converted Hi mass function following the procedure described in the 
text. 



For example, Hoffman et al. [0 conducted a series of pointed observations with the Arecibo 
telescope of known galaxies in the Virgo cluster. They constructed a preliminary HiMF and 
concluded that there is no excess of gas rich dwarf galaxies. Briggs & Rao [^] reanalyzed these 
data supplemented with the Tully & Fisher catalog of H I observations of spiral galaxies, and 
were able to construct the HiMF over the range 10"^ to 10^° M©. They also concluded that there 
is no evidence for a sharp rise in the number of gas rich dwarf galaxies, and that in clusters the 
HiMF might even go down. 

Although the estimation of the space density of gas rich objects using galaxies that are 
selected on their star light is very instructive, it does not yield an unbiased HiMF for all gas 
containing objects. The HiMF determined by optically selected galaxies defines a lower limit 
to the "real" HiMF, and it is not inconceivable that large populations of extragalactic objects 
are missing from the analysis. It has become clear of late that low surface brightness (LSB) 
galaxies contribute significantly to the number density of galaxies in the local Universe (e.g., 
3?, 0). Since these galaxies are found to be rich in HI they are likely candidates to make 



a contribution to the HiMF, while they are easily missed in optical surveys Objects that 
are most certainly missed in surveys in optical wavelengths are hypothetical free floating Hi 
clouds. A conversion from a LF to a HiMF would never include these objects. 



2.2 HiMF Based on Hi Selected Galaxies 



The most ideal method to measure an HiMF is to use a "blind" survey in Hi. No biases 
against galaxies' surface brightnesses or magnitudes are introduced that way. The HiMF can 
be constructed directly and suffers no bias against LSB galaxies, gas rich dwarfs or even free 
floating Hi clouds. The major problem is that finding a galaxy in Hi is much more difficult 
than finding it optically. Blind surveys in H I take hundreds of hours of observing time to yield 
only few dozen galaxies, while optical surveys systematically produce catalogs of thousands of 
galaxies. For illustration, the first published blind Hi survey in the field (Shostak [Q) took 



many days of observing time with the NRAO 300 ft telescope and yielded only one detection 
which later turned out to be a high velocity cloud gravitationally bound to the Milky Way 
galaxy. In order to increase the detection efficiency many surveys were pointed toward known 
overdensities: groups and clusters of galaxies. The M81, Sculptor, CVn I and NGC 1023 groups 
12, 20, 11 as well as the the Hydra [25|, Hercules [0], Centaurus and Fornax |]l| clusters have 
been surveys extensively in the 21cm line. The conclusion that can be drawn from these surveys 
are 1) the shape of the HiMF is similar to that that is derived in the field from optically selected 
galaxies, 2) no excess of gas rich dwarfs has been found, and 3) no Hi clouds without stars 
are found. Also global underdensities (voids) have been surveyed ||39|, 



131, but no Hi selected 



galaxies have been found there. 

A fair comparison of the HiMF for Hi selected galaxies with that from optically selected 
galaxies requires a blind H I survey of the field, with no preference to known over- or underden- 
sities. Henning [|l^, |14| conducted a series of pointings on lines of constant declination over a 
redshift range of 7200 km/s. A total number of 39 significant detections were recorded, of which 
50% were previously unknown. While Hennings HiMF seems to be indicative of an increasing 
number of dwarf galaxies, the overall function lies below the lower limit to the HiMF set by 
counting the optically selected population. 

Two large surveys in the 21cm line have been conducted recently, both with the Arecibo 
telescope. The results from one of these surveys, named AHiSS (Arecibo H l Strip Survey) is 
the topic of the next section. The other survey, similar in size, is discussed in Spitzak 
and Schneider |^T|. There are currently surveys in progress at Dwingeloo [|1^], at Arecibo (The 
Dual-Beam Survey l3l|), and at Parkes 0]. The latter will cover the whole southern sky out 
to 12,700 km/s. 



3 The Arecibo H I Strip Survey 

The survey was carried out during the period of August 1993 until February 1994 when the 
Arecibo telescope was being upgraded. Since the pointing was immobilized, the data were taken 
in drift-scan mode. Two strips of constant declination [6 = 14° and 6 = 23°) were traced, 
together covering 17 hours of RA. The total sky coverage, including the side lobes, was 65 
square degrees, wile the survey depth was cz = 7400 kms~^. The minimal detectable Hi mass 
at the full depth of the survey was 1.5 x 10^ h^QQ Mq in the main beam, while 6 x 10^ /ifoo 
of Hi could be detected at 7 /ifgo Mpc. 

A total of 66 significant (5cr) detections were found, of which 32 are listed in existing catalogs 
of optically selected galaxies. All detections that are sufficiently far away from the Galactic 
equator (6 > 10°) to avoid severe extinction by Galactic dust were imaged in the 5-band with 
the 2.5 Isaac Newton Telescope on La Palma. All Hi signals were identified in the optical 
images, and were found to be ordinary galaxies, having both stars and gas. No free fioating 
H I clouds were detected in this survey, which puts serious constraints on the existence of these 
objects which are potentially a population of protogalaxies. 

21cm follow-up with the VLA in D-configuration was performed on all significant detections. 
These observations were required to determine more accurate fiuxes and positions for all sources. 
Another important reason for doing follow-up with higher angular resolution, is to see whether 
the signals were caused by single clouds or galaxies, or by pairs or small groups, whose line 
emission might stack up in the same channels. We found that this latter situation to occurs in 
five of the 61 cases. For surveys which make use of dishes smaller than the Arecibo telescope, 
and therefore scan the sky with larger beam sizes, this effect is probably more important. Higher 
resolution 21cm follow-up is therefore essential to accurately determine the number of galaxies 



in the survey, their H I masses and to make adequate identifications with optical images. 

From a comparison of the H I and optical properties between the cataloged and uncataloged 
galaxies in our H I selected sample we conclude that the uncataloged galaxies are preferentially 
the ones with lower average optical surface brightness, lower luminosity, lower H I mass, smaller 



linear sizes and higher gas to light ratios However, this does not imply that a blind Hi 
survey yields a population of galaxies that has previously gone unnoticed in optical galaxy 
surveys. Simply because in these two strips the sky has been surveyed deeper in H I than in the 
optical, the uncataloged galaxies will mostly be dwarf galaxies or galaxies with lower optical 
surface brightness. Furthermore, the restricted bandwidth of the receiving system imposes a 
strict limit on the maximum redshift of the detected galaxies, this limits does not exist in 
optical surveys. These two effects cause a predominance of low luminosity systems in an Hi 
selected galaxy sample. It is well established that decreasing luminosity and decreasing surface 



brightness correlate well with increasing Mm/L ratio ^. Therefore, the differences found 
between cataloged and uncataloged galaxies are a natural result of the survey technique and do 
not imply that the new Hi selected galaxies have properties that set them apart from known, 
optically selected galaxies. 

A HiMF for this sample of Hi selected galaxies is determined by using the 1/Vmax method 



IU| . This method consists of summing the reciprocals of the volumes corresponding to the 
maximum distance at which each object could be placed and still remain within the sample. 
Advantages of this method are that the HiMF is automatically normalized and that it is non- 
parametric, that is it does not use a Schechter function as an intrinsic assumption about the 
shape of the HiMF. A recent overview of the different galaxy luminosity function estimators 



is given by Willmer who tests the validity of different methods by means of Monte-Carlo 
simulations. Careful examination of his tables shows that the 1/Vmax method (with binning 
in magnitudes) recovers the input luminosity function satisfactorily, and equally well as the 
more conventional parameterized maximum likelihood method or the stepwise maximum 
likelihood method (SWLM §). 

Fig. shows the HiMF for the AHlSS. The solid points show the HiMF per half decade, 
normalized such that they represent the number of galaxies per Mpc^ per decade. Errorbars 
represent la errors from Poisson statistics. The solid line indicates a best fit Schechter function 
to the points, parameterized by 

e(M„,)«„.r(||)%.p-(^|).(^), ,3) 

where a is the slope of the faint tail, 6* is the normalization, and Mj^j is the characteristic mass 
that defines the knee in the HiMF. Because there are no H I masses lower than 10^ or higher than 
10^° Mq detected by the AHiSS, the HiMF is almost unconstrained in these regions. However, 
upper limits to the space density of H I rich galaxies or H I clouds can be determined by showing 
the sensitivity function in the same figure. The sensitivity function is defined by = 1/V, where 
V is the effective search volume. The arrows indicate upper-limits to the HiMF, on the left they 
follow this sensitivity function, on the right they are determined by a complementary survey 
with Arecibo over the redshift range cz = 19,000 to 28,000 kms~^. The best fit parameters of 
the Schechter function are a = -1.20, 9* = 0.014 Mpc'^ and log(Mj^i/M0) = 9.55. 

3.1 The integral Hi density 

As was mentioned in section |l], the local HiMF can be used to determine the integral gas density, 
fig, at z = 0. At high redshifts Qg can be derived from the study of damped Lja systems, high 
column density gas clouds seen in absorption in the spectra of background quasars. At low 
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Figure 2: H I mass function from the Arecibo H I Strip Survey. The thin hne is the sensitivity 
of the survey defined by = 1/V, where V is the effective search volume. In the region 
IO^Mq < Mhi < 1O^°M0 this function defines an upper limit to the space density of intergalactic 
Hi clouds without stars. The measured Hi mass function per decade is shown by the points. 
The thick line is a Schechter luminosity function with parameters as given in the upper right 
corner. The arrows show upper limits to the space density of H I rich galaxies or intergalactic 
Hi clouds. The two arrows on the right are from a complementary survey with the Arecibo 
telescope over the range 19,000 to 28,000 kms~^. 



redshifts, these measurements are hampered by poor statistics for several reasons: due to the 
expansion of the Universe the expected number of absorbers decreases with decreasing redshift, 
the Lya line is not observable from the ground for redshifts smaller than 1.65, and starlight and 
dust in the absorbing foreground galaxies hinder the identification of the background quasars. 
The most recent determinations of the gas density from damped Lya systems show that the 
redshift range from to 1.65, which corresponds to more than three quarters of the age of the 
Universe (go = 0.5), is covered by only about 20 systems [HO]. Furthermore, the lowest redshift 



damped Lya system currently known is at z = 0.09 illustrating that a z = point can not 
yet reliably be derived from damped Lya observations. 

The gas density detected in the AHiSS can be readily determined by taking the integral over 
the HiMF multiplied by Mhi. This yields pni = r(2 + a)9*M^j^, where pni is the neutral gas 
density in Mq/Mpc^. Using the values from Fig. |^ we derive that pni = 5.8 x 10^ /iioo Mq/Mpc^, 
corresponding to f2g = (2.7 ± 0.5) x 10~^ ^roo- This value corresponds very well with earlier 
estimates based on optically selected galaxies by Rao & Briggs who found pni = 4.8 x 



lO"^ hioo Mq/Mpc . This value has been confirmed by Natarajan & Pettini who have 



used LFs from recent large redshifts surveys. Also the first results of the HiPASS survey [|T7 
yield a similar value. When using their Schechter values a = —1.35, 6* = 0.014 Mpc~^ and 
log(Mj^i/M0) = 9.50, we derive pm = 6.1 x 10^ /iioo Mg/Mpcl 

Figure ^ shows how this gas density is distributed over different H I masses and different lu- 
minosities. As in Fig. |l] the lines represent converted LFs ||23| for different morphological types. 
It is clear that the gas density in the local Universe is dominated by massive galaxies in the 




Figure 3: Left panel: Space density of Hi mass per decade contained in objects of different 
mass. Right panel: Space density of Hi mass per mag contained in objects of different absolute 
magnitude. Lines are as in Fig. |I]. Most of the Hi at the present epoch is locked up in luminous 
galaxies with high H I masses. Dwarfs seem to contribute very little to the gas density. 

range 10^ M0 < Mhi < 10^° M©, with high luminosities. The presently known dwarf galaxies 
are minor contributors to Qg since galaxies with Hi masses below 10^ are responsible for 
only ~ 20% of the gas density. 

Although the run of f2g with redshift is generally considered to reach a peak at z 3 and 
then declines until z = 0, the most recent (preliminary) results from Turnshek et al. show 
that ^g{z) is consistent with no evolution from z = 3 to z ~ 0.15. Still, the 2; = point, 
evaluated via the 21cm in emission, lies approximately a factor of six lower. If, like the AHiSS 
results suggest, most of the neutral gas indeed resides in large spirals, the identification of 



galaxies responsible for damped Lya absorption would be trivial. Rao et al. show that 
the opposite is true: none of the galaxies in the vicinity of the quasar in which two nearby 
absorbers are found are luminous spirals. Also at higher redshifts there are indications that 
dwarf galaxies contribute significantly to the cross section of damped Lya absorbers. Le Brun 



et al. [|l^ find that most of the optical counterparts of damped Lya systems at high redshift 



are dwarf galaxies. Matteucci et al. |22[ conclude on the basis of chemical evolution models 
that most damped Lya absorbers are in fact dwarf galaxies. 

There are two possible solutions to this z = Q conundrum. Either there has been strong 
evolution at very low redshift, implying that most of the star formation has occurred over the 
last few Gyrs, or estimates of the gas density at z = from the 21cm line are still incomplete. 
At least in the range Mhi > 10^ most authors agree well on the shape and normalization of 
the HiMF and derivations based on optically selected and Hi selected galaxies are consistent. 
If an incompleteness of the HiMF exists it must be in the tail of the HiMF. The question to 
answer is therefore: are there large populations of gas rich dwarf galaxies or Hi clouds with 
typical 21cm fluxes undetectable in present surveys? 

In this light it is interesting to note that are recent claims of a sharply rising LF, typically 
faintward of M = —16 (see [pi|]). The galaxies responsible for this upturn are generally of late 



morphological type, with blue colors and low optical surface brightness. These are the galaxies 



normally regarded as being rich in neutral gas (see also ||3T|). 

Another interesting recent result comes from Blitz et al. who argue that many of the 
high velocity clouds (HVCs) populated around the Milky Way Galaxy are actually Hi clouds 
with masses of a few times 10^ distributed trough the Local Group. Their integrated mass 
may account to as much as 10^^ M©. 



4 The Faint Tail 

Obviously, determining the faint tail of the HiMF deserves special attention. An accurate 
determination below Mhi = 10^ has proven to be an extremely difficult exercise in blind 
surveys. For the Parkes Multibeam survey, the largest Hi survey to date, the limiting depth 



for a galaxy with 10 Mq of Hi is 1 Mpc This estimate assumes that the profile width 



of the galaxy is 200 kms~^, while in reality the profile width for these tiny Hi masses will 
probably be much smaller, implying that the same fiux is spread over fewer channels which 
makes the galaxy more easily detectable. Even if the limiting depth is 2 Mpc, the maximum 
search volume for 10^ Mq Hi masses is V = \^D^ ~ 17Mpc^, where Vt is the solid angle. The 
search volume will therefore be confined to the Local Group environment. 

For a blind H I survey, the survey volume as a function of H l mass can be evaluated in the 
following way. The minimal detectable Hi mass MHi^min fo^' ^ So" detection can be expressed as 
236 X ba/SVD^, where a is the noise in the detection spectra in mJy, AV is the profile width in 
km/s and D is the distance to the source in Mpc. Assuming optimal smoothing of the survey 
spectra, the noise goes as cr oc Ay^^/^, and therefore MHi^min oc -D^AF^/^. The limiting depth, 
the maximum distance at which a source could be placed and still be detected, is then given by 
Dc oc Mhi^^^AV^"^/^. The Hi mass is related to the profile width via the Tully-Fisher relation 
and the relation between H I mass and optical luminosity. The "H I Tully-Fisher relation" goes 
as AV cx Mm^^^ The limiting depth is therefore oc Mri^''^^. Finally, since the total 
search volume is V = ^VlD^, we derive V oc M^ii^^. In contrast to optical surveys. Hi surveys 
are conducted with a limited bandwidth which imposes a hard upper limit to the limiting depth 
for massive galaxies. The surveys volume therefore reaches a maximum for H I masses that can 
be detected at the redshift corresponding to the receiver's frequency limit and does not increase 
for higher masses. The sensitivity curve, the inverse of the survey volume, is therefore a straight 
line with slope —1.25, leveling off to a straight line for Hi masses which can be seen all trough 
the survey volume. The decreasing part is parallel to an HiMF with a slope of —2.25. 

The number of detections in a survey can be estimated by taking the integral over the HiMF 
times the survey volume: 

N= / e(MHl)V(MHl)rfMHI. 
JO 

Graphically, this number can be regarded as the surface enclosed by the HiMF and the sen- 
sitivity curve (see Fig. ^. If a < —2.25 the surface would not be closed and the number of 
detections would be infinite. For the HiPASS survey, the instrumental depth limit set by the 
bandwidth of the receiving system is 127 /ij~Qo Mpc and the covered solid angle is = 6.5. If 
then the Schechter parameters for the HiMF found by Zwaan et al. are used, the number 
of detections for this survey will be A/" = 2516 galaxies, increasing the size of existing samples 
by a factor of 30. The statistics will be dominated by Mj^j galaxies: more than 70% of the 
detected galaxies will have Hi masses in the range 10^ M© < Mhi < 10^° M©. If the faint slope 
of q; = —1.2 extends to the lowest masses, only 4 ± 2 galaxies with Hi masses below 10'' M© 
will be detected. This number rises to 13 if the latest value of a = —1.35 quoted by Kilborn 



[p!7| is used. If the slope is much steeper than currently assumed, for example a = —1.7 below 



Mhi = 10^ M0, in concordance with the estimates of the optical LF, the number of dwarfs will 
be 57, high enough to reliably determine the shape of the HiMF down to Mhi = 10^ Mq. 

Thanks to their smaller dishes, synthesis instruments like the WSRT and the VLA sample 
larger volumes per single pointing than the Parkes telescope. These instruments may be more 
useful for constraining the faint end. If the accessible bandwidth is 20 MHz the sampled volume 
will be V = |Dj^a^x(FWHM/2)^7r = 1.5 Mpc^ per single pointing. In 24 hours Hi masses of 
10'' Mq will be detectable throughout the volume. Adopting the HiMF parameters for field 
galaxies from the AHiSS we estimate that the number of detections will be approximately half 
a galaxy per pointing of 24 hours. The statistics will be dominated by galaxies in the range 
10^ M0 < Mhi < 10* M©, two orders of magnitude lower in Hi mass than in the HiPASS 
survey. If again the HiMF with a = —1.7 for Mhi < 10* Mq is assumed, the number of 
detections in the 6 to 6.5 Mq bin will be 0.05. In order to get a statistically significant result 
many months of observing time will be required. Apparently, determining the faint end slope 
of the HiMF below 10^ Mq in the field is not possible with the current generation of radio 
telescopes. 

A way to overcome the poor statistics and the uncertain volume corrections is to study 
nearby groups of galaxies. Typical overdensities of these groups are a factor 25 of the cosmic 
mean, which improves the statistics significantly. Furthermore, the distances of groups can be 
determined quite accurately using Cepheid measurements or by photometry on the brightest 
resolved stars of only one or a few group members, or by using the POTENT method 0. Studies 
of the HI content of galaxies in different environments have shown that the shape of the HiMF 
for Mhi > 10* is independent of cosmic density. This finding justifies the study of the HiMF 
in groups of galaxies, although it is possible that the faint tail may depend on environment. 
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